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We present the relativistic calculation of the hyperfine structure in pionic and kaonic atoms. 
A perturbation method has been applied to the Klein-Gordon equation to take into account the 
relativistic corrections. The perturbation operator has been obtained via a multipole expansion of 
the nuclear electromagnetic potential. The hyperfine structure of pionic and kaonic atoms provide 
an additional term in the quantum electrodynamics calculation of the energy transition of these 
systems. Such a correction is required for a recent measurement of the pion mass. 
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I. INTRODUCTION 

In the last few years transition energies in pionic [l| 
and kaonic atoms have been measured with an un- 
precedented precision. The spectroscopy of pionic and 
kaonic hydrogen allows to study the strong interaction 
at low energies [1, |3, Q by measuring the energy and 
natural width of the ground level with a precision of few 
meV [1,0, Si- Besides, light pionic atoms can addition- 
ally be used to define new low-energy X-ray standards [9] 
and to evaluate the pion mass using high accuracy X-ray 
spectroscopy [H, [111, [l^, El [• Similar endeavour are in 
progress with kaonic atoms 3| • 

In this paper we present the calculation of the hy- 
perfine structure in pionic and kaonic atoms consider- 
ing the perturbation term due to the interaction be- 
tween the pion or kaon orbital moment with the magnetic 
momentum of the nucleus. Non-relativistic calculations 
for the pionic atom hyperfine structure can be found in 
Ref. [ij, 111, m. Other theoretical predictions for HFS 
including relativistic corrections can be found only for 
spin-i nucleus [13, [3 ■ Contrary to these methods, our 
technique is not restricted to this case and it can be used 
for an arbitrary value of the nucleus spin including auto- 
matically the relativistic effects. In particular, we calcu- 
late the HFS energy splitting for pionic nitrogen, which 
has been used for a recent measurement of t he p ion mass 
aiming at an accuracy of few ppm [l^, [HI [3 [l^ ' ^'^'i 
for kaonic nitrogen that has been proposed for the kaon 
mass measurement 2]. 

This article is organized as follows. In Sec. |TT] we cal- 
culate the first energy correction applying a perturbation 
method for the Klein-Gordon equation. In Sec. lIIII we will 
obtain the perturbation term using the multipole expan- 
sion of the nuclear electromagnetic potential. Section HVl 
is dedicated to the numerical calculations for some pionic 
and kaonic atoms. 



II. CALCULATION OF THE ENERGY 
CORRECTION 

The relativistic dynamic of a spinless particle is de- 
scribed by the Klein-Gordon equation. The electromag- 
netic interaction between a negatively charged spin-0 par- 
ticle with a charge equal to q ^ — e and the nucleus 
can be taken into account introducing the nuclear po- 
tential in the KG equation via the minimal coupling 
Pv ^ Pu — qAi, [20]. In particular, in the case of a cen- 
tral Coulomb potential (Vo(''),0), the KG equation for a 
particle with a mass m is: 



1 



[ihdt + eVo{r)Y + h'V'\M>o{x), (I) 



where h is the Planck constant, c the velocity of the light 
and the scalar wavefunction ^'o (x) depends on the space- 
time coordinate x = {ct, r). We consider here the station- 
ary solution of Eq. ([1]). In this case, we can write: 



"^oix) = cx.-p{-iEot/h) ipair) 
and Eq. ([T|) becomes: 
I 



(2) 



[Eo + eVa{r)Y + rV^ - m'c' } (^o(r) = 0, (3) 



where Eq is the total energy of the system (sum of the 
mass energy mc^ and binding energy So). 

The perturbation correction Ei can be deduced in- 
troducing an additional operator W in the zeroth order 
equation: 
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[E + eVo{r)f + h^V^ 



m^c^ - W(r) } ip{r) = 



(4) 

W{r) is in general non-linear. In the case of a correction 
Vi to the Coulomb potential Vb, we have: 



W(r) 



1 



(2e2yo(r)Fi(r)-|-2eS Vi(r)) (5) 
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If we consider the interaction with the nuclear magnetic 
field as a perturbation, we have: 

W{r) = ihe [2A,{r)d' + [d^,A'{r)]] - e^A'{r)Adr). (6) 
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The correction to the energy due to W can be calcu- 
lated perturbatively with some manipulation of Eqs. ([3]) 
and (UJ [U, [13] , or via a linearization of the KG equation 
using the Feeshbach- Villars formalism [2^ [23| • In both 
cases we have 



El 

where we define 

{A) = . 



c2(VF) 



2 ( K"' 



(eVn 



(7) 



(8) 



Equation ([7]) is valid for any wavefunction normalization. 



III. CALCULATION OF THE HYPERFINE 
STRUCTURE OPERATOR 

The expression for W{y) in the HFS case is derived 
using the multipole development of the vector potential 
A(r) in the Coulomb gauge [H, [H, [l^l • We neglect here 
the effect due to the the spatial distribution of the nuclear 
magnetic moment in the nucleus [ll] (Bohr-Weisskopf 
effect), while effect due to the charge distribution (Bohr- 
Rosenthal effect) are included in the numerical results of 
SecHVl 

The hyperfine structure due to the magnetic dipole 
interaction is obtained by taking into account the first 
magnetic multipole term. Using the Coulomb gauge we 
have in,!!!: 



A(r) = -~i^V2r 
An 



(9) 



where the symbol "o" indicates here the general scalar 
product between tensor operators. operates only on 
the nuclear part |/m/) and C^^ is the vector sperical 
harmonic [26|, acting on the pion part \nlm) of the 
wavefunction. We can decompose the perturbation term 
W{r) as: 



W{v)^Wi{v) + W2{r), 



(10) 



where 



W^i(r) = +ihe [2A,{Y)d' + [9„ A'(r)]] , (11) 

is the linear part and 

W2{y) = -e^A\v)A,{v), (12) 

is the quadratic part. 

We study first the operator Wi. We note that 
[9i,yl*(r)] = —V • A(r) = since we are using the 
Coulomb gauge. In this case we have: 

W^i(r) = +2ineA,{Y)d' = -2i?ieA(r) • V = 

_e^o?i-^r-2(c" . V)oMi. (13) 
27r 



Using the properties of the spherical tensor [26|, l29| , we 
can show that: 



— L 

V2 



(14) 



where Lq is the dimensionless angular momentum oper- 
ator in spherical coordinates. The perturbation operator 
can be written as a scalar product in spherical coordinate 
of the operator acting on the pion wavefunction, and 
the nuclear operator M^: 



1^1 (r) 



e^Qh 
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'(ioM^) = oM^ (15) 



with 



1 _ eMO?i 3r 



(16) 



The expected value of the operator Wi can be evalu- 
ated applying the scalar product properties in spherical 
coordinates f29l. [30|: 



{n'l'\\T^nl){I\\M^\\I), (17) 



{n'l'lF'm'p\Wi\nlIFmF) = (-1)'+'+^ (5ff'^™^„,'^<5//' X 

F I I' 
1 / / 

where { ^ g j } represents a Wigner 6-j symbol. The 

reduced operator {n'l'\\T^\\nl) is calculated from the ma- 
trix elements {n'l'm'\T^\nlm) by a particular choice of 
the quantum numbers m and q applying the Wigner- 
Eckart theorem: 



(n'l'WT'Wnl) = Sio 



(-1) 



i-i 



I 1 I 
-1 1 



{n'l'l\T^\nll) = 



SioVlVTTTV2rTT^^{nl'l\r-^L,\nll) 



2n 



= SioSwVlVTTlV2rTl^(n'l\r-'\nl), (18) 

ZTT 

where ( ^ g j ) indicates the Wigner 3-j symbol. 

The nuclear operator can be related to the magnetic 
moment on the nucleus by {II\Mq\II) = fijfiN [29. [27| 
where /z/ is the nuclear dipole momentum in units of the 
nuclear magneton /i^v — eh/2mpc: 



{I\\M'\\I) = 



I 1 I 

-I I 



(19) 



Considering Eq. the total expression for Tyi(r) 

becomes: 

{nl'lF'm'p\Wi\nlIFmF) = 
eiioh F{F + 1) - /(/ + 1) - /(/ + 1)) 



27r 



21 



{n'l\r~^\nl). (20) 
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which, as expected, is equal to zero for I — (then / — 
F). 

To find the final expression of the HFS energy shift, we 
have to evaluate the contribution of the operator W2 (r) = 
—e^A^{r)Ai{r) in the (W) diagonal terms. Using Eq. ([9]), 
we have: 

(nlIFmF\W2\nlIFmF) = +2 (y^)^ x 

{nlIFmp\{r-'^C^^ o M^) • (r^^c" o M^)\nlIFmF). 

(21) 

We are in presence of three independent scalar products: 
two scalar products between the tensor C^^ and the vec- 
tor M^, and the scalar product between the vectorial op- 
erators C^^ o M^. The "•" scalar product in W2 can be 
decomposed using the properties of the reduced matrix 
elements of a generic operator product , of rank K, 
between non-commutating tensor operators and V'' 
of rank k. For our case, this scalar product corresponds 
to a tensor product with K = 0, and fc = 1: 

X° ^U^-V^ = (r-^c" o Ml) • (r-^c" o M^), (22) 
f/i = yi = (r^^C" o M^). (23) 

We have fH: 



Applying the Wigner-Eckart theorem we obtain 

{nlIF'F'\V^\nlIFF) = 

[nl\\r-^C^\\nl){I\\M^\\I)x 



r2V2 



^ {lm'Im'j\lIF'F'){lniInii\lIFF): 



I 1 I 

-to' 1 m 



I 1 I 

-to' —1 TO 



I 1 / 

—m'j —1 TO/ 



1 



-1 TO/ 



(28) 



The reduced matrix element {nl\\r '^C^\\nl) can be de- 
composed in a radial and angular part 



{nl\\r-^C^\\nl) = (n/|r-2|„^)(;||ci||;). 



(29) 



Due to the symmetry properties, (/||Ci||Z) is equal to zero 
for any I [2^ . This result implies that the reduced matrix 
elements of and are always equal to zero. As a 
consequence, the diagonal elements {A^{r)Ai{r)) = for 
any wavefunction, i.e., W2 does not contribute to the 
HFS energy shift. 



{nlIF\\X^\\nlIF) =^ 



1 1 

F F F' 



{nlIF\\U^\\nlIF') {nlIF'\\V^ \ \nlIF). (24) 

V' and are scalar products between commuta- 
tive tensor operators, and their reduced matrix element 
(F'||F||i^) can be calculated applying again the Wigner- 
Eckart theorem for the component q = Q: 

-F' F 

where |3l|: 



We can now write the final expression for the HFS 
energy correction: 



E' 



nlF 



/i//iAre/io?ic2 



{nl\Vo{r)\nl)) 
F{F + 1)-I{I + 1)- 1(1 + 1] 
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{nl\r-'^\nl). (30) 



This formula is obtained by a perturbation approach 
of the KG equation. For this reason, all the relativistic 
effects are automatically included in Eq. (I5D1) . In the 
non-relativistic limit c — > cx), {Eq — {V))/c'^ — > to and we 
find the usual expression of the HFS for the Schodinger 
equation [s^l 



1/0' -C/o' -^~'(C"oMi)o 



V2 ' 



(26) 

To evaluate the matrix element (nllF' F'\V(^\nlIFF) 
we can explicitly decompose \nlIFF) as a function of 
the eigenfunctions \nlm) and |/to/) using the Clebsch- 
Gordan coefficients. We have 

(nlIF'F'\Vo^\nlIFF) 

— ^ V {lm'Im'j\lIF'F'){lmImi\lIFF)x 

[{nlm'\r-^Cl\nlm){Im'j\M^j^\Imi)- 

{nlm'\r~^C\\nlm){Im'j\Ml\Imi)] . (27) 



IV. NUMERICAL RESULTS AND BEHAVIORS 

We present here some calculations for a selection of 
pionic and kaonic atom transitions. Such calculations are 
obtained solving numerically the Klein-Gordon equation 
using the multi-configuration Dirac-Fock code developed 
by one of the author (P.I.) and J. -P. Desclaux [33l. |34|. 
Issl . [3^ that has been modified to include spin-0 particles 
case, even in the presence of electrons [131 . The first part 
is dedicated to the 5^4 and 8 — > 7 transitions in pionic 
and kaonic nitrogen, respectively. In the second part we 
will study the dependence of the HFS splitting against 
the nuclear charge Z to observe the role of the relativistic 
corrections. 
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TABLE L Energy (in eV) contribution for the selected levels 
in pionic nitrogen. The first error takes into account neglected 
next order QED corrections. The second is due to the accu- 
racy of the pion mass (±2.5 ppm). 





5g-4f 


5f-4d 


Coulomb 


4054.1180 


4054.7189 


Finite size 


0.0000 


0.0000 


Self Energy 


-0.0001 


-0.0003 


Vac. Pol. (Uehling) 


1.2485 


2.9470 


Vac. Pol. Wichman-KroU 


-0.0007 


-0.0010 


Vac. Pol. Loop after Loop 


0.0008 


0.0038 


Vac. Pol. Kallen-Sabry 


0.0116 


0.0225 


Relativistic Recoil 


0.0028 


0.0028 


HPS Shift 


-0.0008 


-0.0023 


Total 


4055.3801 


4057.6914 


Error 


±0.0011 


±0.0011 


Error due to the pion mass 


±0.010 


±0.010 



A. Calculation of the energy levels of pionic and 
kaonic nitrogen 

The precise measurement of 5g 4/ transition in pi- 
onic nitrogen and the related QED predictions allow for 
the precise measurement of the pion mass [Tl|, [H, [H, 
[l^. In the same way, the transition 8k 7i in kaonic 
nitrogen can be used for a precise mass measurement of 
the kaon ^2]. For these transitions, strong interaction ef- 
fects between meson and nucleus are negligible, and the 
level energies are directly dependent to the reduced mass 
of the atom. The nuclear spin of the nitrogen isotope 
^'^N is equal to one leading to the presence of several 
HFS sublevels. The observed transition is a combination 
of several different transitions between these sublevels, 
causing a shift that has to be taken into account to ex- 
tract the pion mass from the experimental values. Tran- 
sition probabilities between HFS sublevels can easily be 
calculated using the non-relativistic formula [s^, HI] (the 
role of the relativistic effects is here negligible), if one 
neglect the HFS contribution to the transition energy: 



lIF^n' I'lF' 



(2F+1)(2F' + 1) \ V F' I 



2/+ 1 



F 



F' I) 

I ^ j Anl^n'l', (31) 



where 



_ 4(g„; - g„> ;_i)3 a I ,,,,, 

^nl-*n'l' — o 9 /H- '^(;)2 2Z±1 ' 



Sm^c'^h {Za) 



with 



(32) 



Tyn'V _ 
^nl — 



^0 "'0 



Ki{r)^n'i'{r)r''dr. (33) 



Co = h/{mcZa) is the Bohr radius and are the non- 
relativistic wavefunctions. 

For these calculations, presented in Tables HI arid HH we 
used the nitrogen nuclear mass value from Ref. [39| . The 



TABLE II: Hyperfine transition energies and transition rate 
in pionic nitrogen. 

Transition F-F' Trans, rate (s"^) Trans. E (eV) Shift (eV) 



5^ _> 4^; 


4-3 


4.57 




10" 


4057 6876 


-fl 00606 




3-2 


3.16 


X 


10" 


4057.6970 


0.00341 




3-3 


2.98 


X 


10" 


4057.6845 


-0.00910 




2-1 


2.13 


X 


10" 


4057.7031 


0.00946 




2-2 


2.25 


X 


10" 


4057.6948 


0.00112 




2-3 


0.01 


X 


10" 


4057.6822 


-0.01138 


5g 4d 


5-4 


7.13 


X 


10" 


4055.3779 


-0.00304 




4-3 


5.47 


X 


10" 


4055.3821 


0.00113 




4-4 


5.27 


X 


10" 


4055.3762 


-0.00482 




3-2 


4.17 


X 


10" 


4055.3852 


0.00420 




3-3 


0.36 


X 


10" 


4055.3807 


-0.00029 




3-4 


0.01 


X 


10" 


4055.3747 


-0.00624 



Coulomb term in the Table includes the non-relativistic 
recoil correction using the reduced mass on the KG equa- 
tion. The pion and nucleus charge distribution contribu- 
tion are also included. The pion charge distribution ra- 
dius contribution is included following [37|, S3| ■ For the 
pion charge distribution radius we take = 0.672 ±0.08 
|4ll |. For the nuclei we take values from Ref. The 
leading QED corrections, vacuum polarization, contribu- 
tion is calculated self-consistently, thus taking into ac- 
count the loop-after-loop contribution to all orders, at 
the Uehling approximation. This is obtained by including 
the Uehling potential into the KG equation [3a]. Other 
Higher-order vacuum polarization contribution are cal- 
culated as perturbation to the KG equations: Wichman- 
KroU and Kallen-Sabry [H, |43| . The self-energy is cal- 
culated using the expression in Ref. [i^ and it includes 
the recoil correction. The Relativistic recoil term has 
been evaluated adapting the formulas from Refs. [H, l46j 
(more details can be found in Ref. The calcula- 

tions presented here do not take into account second or- 
der recoil effects (Fig. [T] top), or higher QED corrections 
as vacuum polarization and self-energy mixed diagrams 
(Fig. [H botton). The contribution from these terms has 
been estimated using the formula for a spin-i particle 
with a mass equal to the pion's. For the 5 — > 4 pionic 
nitrogen transitions, vacuum polarization and self-energy 
mixed diagrams contribute in the order of 1 meV for the 
diagram with the vacuum polarization loop in the nuclear 
photon Hue (Z^ (Fig. [T] bottom-left), and 0.0006 mcV for 
the diagram with the vacuum polarization loop inside the 
self-energy loop [4^ (Fig. [T] bottom- right). The second 
order recoil contributions are in the order of 0.04 meV 
^49*1 (Fig. [T]top). The largest contribution comes from 
the unevaluated diagram with the vacuum polarization 
loop in the nuclear photon line 

Assuming a statistical population distribution of the 
HFS sublevels, we can use Eq. (|31|1 to calculate the 
mean value of the transitions using the results in Ta- 
ble [ni Comparing this calculation with the one without 
the HFS, we obtain a value for the HFS shift. For tran- 
sitions 5.g — > 4/ and 5/ Ad we obtain shifts of 0.8 
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TABLE IV: Hyperfine transition energies and transition rate 
in kaonic nitrogen. 

Transition F-F' Trans, rate (s'^) Trans. E (eV) Shift (eV) 
8i -> 7h 



7i 



7-6 


I 


.19 




10" 


2970 


.4169 


-0, 


.00216 


6-5 


1, 


.00 


X 


10" 


2970, 


.4196 


0, 


.00050 


6-6 


0, 


.98 


X 


10" 


2970, 


.4145 


-0, 


.00453 


5-4 


0, 


.84 


X 


10" 


2970, 


.4217 


0, 


.00265 


5-5 


0, 


.03 


X 


10" 


2970, 


.4175 


-0, 


.00154 


5-6 


0, 


.00 


X 


10" 


2970, 


.4125 


-0, 


.00656 


8-7 


1, 


.54 


X 


10" 


2969, 


.6365 


-0, 


.00149 


7-6 


1, 


.33 


X 


10" 


2969, 


.6383 


0, 


.00029 


7-7 


1, 


.31 


X 


10" 


2969, 


.6347 


-0, 


.00326 


6-5 


1, 


.15 


X 


10" 


2969, 


.6398 


0, 


.00178 


6-6 


0, 


.03 


X 


10" 


2969, 


.6367 


-0, 


.00126 


6-7 


0, 


.00 


X 


10" 


2969, 


.6332 


-0, 


.00480 



FIG. 1: Diagrams relative to the unevaluated QED contri- 
butions: second order recoil correction (top), and vacuum 
polarization and self-energy mixed diagrams (bottom). Their 
effects are estimated using the available formulas for spin-i 
particles. 



act cancellation between the weighted excited sublevels 
energy shifts as seen from Eq. ([31]) . 



B. 



General behavior of the hyperfine structure 
correction over Z 



TABLE III: Energy (in eV) contribution for the selected lev- 
els in kaonic nitrogen. The first error takes into account ne- 
glected next order QED corrections. The second is due to the 
accuracy of the kaon mass (±32 ppm). 





8k- 7i 


8i-7h 


Coulomb 


2968.4565 


2968.5237 


Finite size 


0.0000 


0.0000 


Self Energy 


0.0000 


0.0000 


Vac. Pol. (Uehling) 


1.1678 


1.8769 


Vac. Pol. Wichman-KroU 


-0.0007 


-0.0008 


Vac. Pol. Loop after Loop 


0.0007 


0.0016 


Vac. Pol. Kallen-Sabry 


0.0111 


0.0152 


Relativistic Recoil 


0.0025 


0.0025 


HFS Shift 


-0.0006 


-0.0008 


Total 


2969.6374 


2970.4182 


Error 


0.0005 


0.0005 


Error due to the kaon mass 


0.096 


0.096 



and 2.2 meV, respectively. These values correspond to a 
correction to the pion mass between 0.2 and 0.6 ppm. 

The transition energies for the 8^7 transitions in 
kaonic nitrogen are presented in Tables IIIII and IIVI As 
for the pionic nitrogen, the error contribution due to the 
QED correction not considered is dominated by the un- 
evaluated diagram with the vacuum polarization loop in 
the nuclear photon line [13], the associated correction is 
estimated in the order of 0.5 meV. For 8fc — > 7i and 
Si 7h transitions we have a HFS shift of 0.6 and 
0.8 meV, respectively, which correspond to a correction 
of the kaon mass between 0.2 and 0.3 ppm. 

As a general note, we remark that if we assume a sta- 
tistical distribution of the initial state sublevels popula- 
tions, transitions nl — )■ n' s with a s orbital as final state 
have an average HFS shift equal to zero due to an ex- 



For the non-relativistic case, the HFS splitting normal- 
ized to the binding energy and to the nuclear magnetic 
moment, depends linearly on Za. Any deviation from 
this linear dependence in the Klein-Gordon HFS can be 
attributed only to relativistic effects. 

To study the behavior of the normalized HFS splitting 
(-^^=3/2 ~ ^fL3/2)/(^o/^^) ^'^^ ^^'^ relativistic case, we 
calculated the HFS for a selected choice of pionic atoms 
with a stable nucleus of spin 1/2. The orbital 9p has 
been chosen to minimize the effect of the finite nuclear 
size and strong interaction shifts, particularly for high 
values of Z. The results are summarized in Table IVl For 
these calculations we used the nuclear mass values from 
Ref. [3^, the nuclear radii from Refs. [43, [SOi] and the 
nuclear magnetic moments from Ref. |5l| . 

For higher Z values a non- linear dependence on Za ap- 
pears as we can see in Fig. [51 This non-linearity originates 
in the two different parts of Eq. pO)) : the non-trivial de- 
pendency on Sq in the denominator and the expectation 
value {nl\r^^\nl). 



V. CONCLUSIONS 

We presented a relativistic calculation of the hyperfine 
structure in pionic and kaonic atoms. The precise evalu- 
ation of the specific case of pionic and kaonic nitrogen is 
particularly important for the new measurement of the 
pion and kaon mass. The small error on the theoretical 
predictions, of the order of 1 meV for the 5^4 transi- 
tion, corresponds to a systematic error of ^ 0.2 ppm for 
the pion mass evaluation, considerably smaller than the 
error of previous theoretical predictions [53 |. 
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TABLE V: HPS separation of the F=l/2 and F=3/2 levels 
for the 9p orbital for pionic atoms with spin ^ nucleus. 

Element Z 9p energy (eV) HPS splitting (cV) 



XT 

n 


1 
± 




U.UUUl 


xie 




z 


1 7/1 QQ7n 

-1 / 4.00 ( U 


-u.uuuy 


13(3 


6 


-1633.402 


0.0060 




7 


-2226.813 


-0.0039 


19p 


9 


-3689.435 


0.0767 


31 p 


15 


-10286.63 


0.1544 


5^Fe 


26 


-31023.63 


0.0643 


^^Se 


34 


-53146.21 


0.8293 


89 Y 


39 


-69987.00 


-0.3143 




47 


-101681.2 


-0.4266 


i29Xe 


54 


-134178.0 


-4.1295 




74 


-250634.4 


1.2629 


202pj^ 


82 


-306731.6 


7.8662 



The formalism presented in this article can be applied 
for other effect as the quadrupole nuclear moment which 
can not be negligible for mesonic atoms with high Z . In 
this case, HFS due to the quadrupole moment can be 
predicted using the next multipole in the development 
of the electric potential of the nucleus to evaluate the 
correspondent perturbation operator. This application is 
particularly important for the calculation of the atomic 
levels in heavy pionic ions, where relativistic and nucleus 
deformation effects can be taken into account at the same 
time. 




30 40 50 60 
Atomic number Z 

FIG. 2; Value of relative splitting (S^'Lg/j ~ Ef'=3/2)/(^oI^i) 
for pionic atoms with different values of Z. Nuclear masses, 
radii and ma gne tic moments values have been obtained from 
Refs. [3i,|42,M[5l[. 
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